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Abstract. Heavy quark structure functions from HERA provide a direct handle on the medium and 
small-x gluon PDF. In this contribution, we discuss ongoing progress on the implementation of the 
FONLL General-Mass scheme with running heavy quark masses, and of its benchmarking with the 
HOPPET and OpenQCDrad codes, and then present the impact of the recently released combined 
HERA charm production cross sections in the NNPDF2.3 analysis. We find that the combined charm 
data contribute to constraining the gluon and quarks at small values of Bjorken-x. 
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Charm structure function data and PDF fits. Charm production in deep-inelastic 
scattering is directly sensitive to the gluon PDF. The ZEUS and HI collaborations at 
the HERA collider have measured charm production in DIS with a wide variety of 
techniques, and these F 2 C data is included in all modern PDF fits (see [1] for a recent 
overview). The HERA experiments have recently released their combined data on charm 
production cross sections from HERA Runs I and II [2], where a common consistent 
data set with the full correlation matrix is provided, and in addition the combination 
procedure yields systematic errors rather smaller than what one would expect from the 
naive combination of all the data, because of the mutual cross-calibration between HI 
and ZEUS. 

In this contribution we explore the impact of replacing the separated HI and ZEUS 
F£ data with the combined charm production cross sections o^ c in the NNPDF2.3 
analysis [3]. We also discuss how one can generalize the FONLL General-Mass VFN 
scheme to include running heavy quark masses in the MS scheme. Heavy quark structure 
functions in the MS scheme lead to a better behaved perturbative expansion than pole 
masses [4] and allow to compare the value of the heavy quark masses used with those 
determined by other experiments (e. g. by LEP data). 

FONLL Structure functions in the MS scheme. Treating heavy quark structure func- 
tions in the MS scheme results in an improved convergence of the perturbative expansion 
and allows to consistently compare the mass of the charm quark used in the PDF analysis 
with other determinations. The NNPDF fits are based on the FONLL GM-VFN scheme 
for heavy quark structure functions [5], with pole masses as default. It can be shown that 
FONLL can be extended to use MS heavy quark masses, we discuss here some progress 



TABLE 1. Relative differences at Q 2 = 10 4 GeV 2 for NNLO PDF evo- 
lution in the massless scheme with MS running heavy quark masses as 
implemented in the FastKernel framework in comparison to HOPPET, 

£rei = | ($P(*, Q 2 ) ~ q)° P {x, 2 2 )) /qf{x, Q 2 )\, for various PDF flavor com- 
binations. The PDFs and the settings of the comparison are the same as in 
the Les Houches benchmark comparisons [7]. 
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in the this direction. Let us recall that FONLL can combine different perturbative orders 
for massive and massless structure functions: FONLL-A combines the NLO massless 
with ff(a s ) massive pieces, FONLL-C achieves the same at G (oc s 2 ) and FONLL-B is 
an intermediate case where the & (a 2 ) massive calculation is combined with the NLO 
massless result. 

The first ingredient to consider, as compared to pole masses, is the NNLO massless 
PDF evolution, which is modified by the scheme transformation between pole and MS 
running masses (NLO evolution is identical in the two schemes). In particular, the heavy 
quark thresholds require different matching conditions with running masses as compared 
to pole masses [4]. To benchmark the implementation of the MS NNLO PDF evolution 
in the NNPDF FastKernel framework, we report in Table 1 the relative accuracy 
against the x-space evolution code HOPPET [6], that also has the option to perform 
PDF evolution with running heavy quark masses. This benchmark has been performed 
using the same parameters as in the Les Houches PDF comparison [7]. We find excellent 
agreement over all the kinematic range. 

A second ingredient of the FONLL MS masses implementation is the comparison 
with the OpenQCDrad code 1 , which provides charm structure functions in the Nf = 3 
FFN scheme both for pole masses and for MS running masses. The comparison is done 
as follows. First of all, we compute charm structure functions in the massive scheme 
with pole masses set to M c = \/2 GeV in both codes, and check that there is reasonable 
agreement. Then we transform the pole mass to the MS running mass, which is this 
case corresponds to m c (m c ) = 1.06 GeV, and use this as input of the MS running mass 
computation both in FastKernel in the FFN scheme and in OpenQCDrad. In Table 2 
we present the percentage difference for charm structure function F 2 C in the MS scheme, 
we find good agreement over a wide kinematical range. 

With these two ingredients, it is possible to generalize FONLL to include running 
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TABLE 2. Relative differences between the FastKernel 
implementation of the charm structure function F 2 C with MS 
heavy quark running masses and the OpenQCDrad results, 
erei = \ {F^(x,Q 2 ) - F°^(x,Q 2 )) /F^(x,Q 2 )\. The computation has 
been performed in the massive scheme at order G (a 2 ) with the LH toy 
PDFs. 
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FIGURE 1. The ratio of the gluon (left plot) and quark singlet (right plot) between the NNPDF2.3 
PDFs and the fit where with the separated F% HI and ZEUS data have been replaced by the combined 
(J^ c charm reduced cross sections. PDFs have been evaluated at a typical LHC scale of Q 2 = 10 4 GeV 2 . 



heavy quark masses. Of course the difference between pole and running masses appears 
only at (aj) , so FONLL-A will not be changed, while both FONLL-B and FONLL- 
C will be modified by the scheme transformation. A more detailed discussion of the 
implementation of running heavy quark masses in FONLL, and the corresponding 
impact on parton distributions, will be presented elsewhere. 

The impact of the combined HERA charm production data. The NNPDF analy- 
sis have included all the available HI and ZEUS charm structure function data since 
NNPDF2.1 [8, 9, 10, 3]. We have now implemented the combined HERA charm pro- 
duction cross sections <j£c in the NNPDF code, and performed various fits, at NLO and 
NNLO, with the same settings as NNPDF2.3 but replacing the separated HI and ZEUS 
F£ data with the combined HERA data. We take fully into account the 43 sources of 
correlated systematics of the combined dataset, which include normalization and proce- 
dural uncertainties. To ease the comparison, structure functions are computed in the pole 
mass scheme with the same mass values as in NNPDF2.3. We obtain a good description 
of the data, with a X 2 P er data point of about 1.2. 

We show the impact of the combined HERA charm data on NNPDF2.3 NNLO in 
Fig. 1 and we also show the associated distances between the two fits in Fig. 2. As we 
can see, the net effect of replacing the separated F 2 C data with the combined reduced cross 
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FIGURE 2. Distances between the two fits of Fig. 1. See [11] for the definitions. 



sections is to shift the central value of the small-x gluon and total quark singlet PDFs by 
a moderate amount, between one third and half a sigma of the PDF uncertainty. 

Outlook. We have discussed the implementation of the combined HERA charm data 
in the NNPDF framework, showed that these new data provide some useful constrains 
on the poorly known small-x gluons and quarks, and reported on the implementation of 
the FONLL GM scheme with running heavy quark masses. We plan to use these results 
to perform a determination of m c (m c ) from the combined HERA charm production data, 
with the same techniques used to determine a s (M z ) [12, 13]. 
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